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Assessment of hypnotic susceptibility is usually obtained through the application of psychological 
instruments. A satisfying classification obtained through quantitative measures is still missing, although it 
would be very useful for both diagnostic and clinical purposes. Aiming at investigating the relationship 
between the cortical brain activity and the hypnotic susceptibility level, we propose the combined use of two 
methodologies - Recurrence Quantification Analysis and Detrended Fluctuation Analysis - both inherited 
from nonlinear dynamics. Indicators obtained through the application of these techniques to EEG signals of 
individuals in their ordinary state of consciousness allowed us to obtain a clear discrimination between 
subjects with high and low susceptibility to hypnosis. Finally a neural network approach was used to perform 
classification analysis. 

Classifying individuals according to their level of susceptibility to hypnosis has become necessary within the 
perspective of experimental hypnosis in order to predict individual responses to suggestions. The differ- 
ences in susceptibility are particularly important because they influence individuals' everyday life in their 
ordinary state of consciousness. In fact, out of hypnosis, individuals can also differ in various aspects of sensor- 
imotor integration 14 and cardiovascular control 5 . 

Several scales 6 can be used to characterize subjects as high (highs), medium (mediums) and low (lows) hypno- 
tizables. Nonetheless, the administration of scales requires time and experience; moreover, most of the hypnotiz- 
ability scales measure suggestibility after hypnotic induction rather than hypnotizability, i.e. suggestibility in the 
absence of hypnotic induction 7 . Suggestibility while awake may change after hypnotic induction 8 , and its assess- 
ment is important because hypnotizability predicts the individuals' response to suggestions out of hypnosis and 
allows, for instance, to ease manageability of suggestions such as analgesia and anaesthesia in clinical contexts. 

The few attempts pursued to extract information on susceptibility levels from electroencephalograms (EEGs) 
were performed after hypnotic induction. Many techniques inherited from nonlinear time series analysis and 
nonlinear dynamics have been applied to EEG signals 9,10 , proving the existence of a significant relationship 
between hypnotizability and a few indices extracted from the EEG. At present, and to the best of our knowledge, 
only one study proves that it is possible to classify highs and lows using EEGs recorded in the ordinary state of 
consciousness during relaxation 11 . The difference between highs and lows, however, was present only in the 
earliest minutes of relaxation, in line with earlier reports indicating that highs and lows process the relaxation 
request through different cognitive strategies leading to the same perceived relaxation 12 . In Madeo etal. 11 the EEG 
dynamics of individuals with different levels of hypnotizability was characterized through the quantitative 
measures of two indicators called Determinism and Entropy - related to the determinism and the complexity 
of the EEG signal, respectively - in the framework of the nonlinear time series analysis based on recurrences 13 . 
Actually Determinism and Entropy were recognized as the most appropriate parameters for discriminating highs 
and lows subjects based on a minimal set of quantitative indicators 11 . 

Presently our objective is more refined, aiming at developing an optimal classification tool that can be used for 
clinical purposes, which should be as simple and user-friendly as possible and to avoid the use of a large number of 
indicators that could increase the complexity of the classification algorithm. As a consequence, we propose the 
combined use of Recurrence Quantification Analysis (RQA) and Detrended Fluctuation Analysis (DFA). 
Specifically, in order to obtain a more detailed and broad-spectrum investigation of the EEG recordings, we 
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integrate determinism with the fluctuation exponent as a comple- 
mentary indicator related to the stochastic part of the signal. 

The RQA is a nonlinear technique that can be traced back to the 
work by Poincare 14 , and Ruelle and Takens 1516 . It quantifies the 
small-scale structures of recurrence plots which present number 
and duration of recurrences of a dynamical system in a reconstructed 
phase space. The main advantage of this kind of analysis is to provide 
useful information even for non-stationary data where other meth- 
ods fail. The DFA is a method basically designed to investigate long- 
range correlations in nonstationary series 1719 , through the estimate 
of a scaling exponent obtained from the slope of the so-called fluc- 
tuation function F(s) as a function of lag s, in a log-log plot. The value 
of this exponent can discriminate between (short or long-range) 
correlated and uncorrelated time series. The observed signals char- 
acterizing a complex system often exhibit long-range correlations. It 
is of crucial importance and significance to quantify such long-range 
correlations in order to have a deep understanding of the dynamics of 
the underlying complex systems. The DFA technique has been exten- 
sively used on EEG databases to investigate the relationship between 
hypnotizability and brain activity 20 , but typically it has been applied 
to hypnotized subjects only and never on datasets obtained from 
basal EEG out of hypnosis. 

As many other physiological phenomena, the hypnotizability of a 
subject emerges from a complex mechanism, involving many cereb- 
ral and cognitive activities, such as memory retrieval, self-reflection, 
mental imagery, and stream-of-consciousness processing 21 . As a 
consequence, the dynamics of the underlying system must be 
explored at many different time scales. This goal can be achieved 
through the integration of the proposed techniques, strengthening 
the results provided by the separate application of the two methodo- 
logies. On one hand RQA is a powerful discriminatory tool which can 
provide information regarding the degree of determinism character- 
izing the system, as well as the degree of complexity or randomness of 
the signals. On the other hand, DFA is used to quantify the fractal- 
like scaling properties of the same signals. Once combined with the 
other indicator estimated in the framework of RQA, the power-law 
exponent of DFA can be used to classify hypnotic susceptibility, 
providing a quantitative discrimination between highs and lows. 

The aim of this study is to show that the use of RQA coupled with 
scaling analysis can discriminate highs from lows on the basis of EEG 
recordings obtained from non-hypnotized participants during 15 
minutes of relaxation 11 . In order to achieve a systematic integration 
of the two methodologies several neural networks were developed 
based on the time varying measures RQA and DFA evaluated on a 
dataset including the EEG recordings of 8 highs and 8 lows. On the 
basis of previous findings 11 four channels - out of 32 - were used in 
the neural network. The neural network approach, unlike typical 
classification or clustering techniques such as linear separation, 
allowed us to take into account the intrinsic nonlinearities present 
in the signals. 

Results 

In this section the preliminary step performed for data acquisition, 
and the main results obtained through the combined use of RQA, 
DFA and Neural Networks, are reported. 

Experimental Setup and Data Preprocessing. The experiments 
were performed on 16 healthy males (age: 19-30 yrs), 8 highs 
(score > 8/12) and 8 lows (score < 3/12), selected according to the 
Stanford Hypnotic Susceptibility Scale 22,23 . Subjects were asked to 
relax their best, remain silent and avoid movements for a time 
period of 15 minutes. The EEG signals have been acquired through 
a Neuroscan device (40 channels) with a sampling frequency of 
1 kHz during the whole experimental session. After acquisition, 
data was preprocessed in order to remove artifacts and measure- 
ment errors. Among the 40 recorded channels, CZ, TP7, P3 and 



POl incorporated the best information to discriminate between 
highs and lows on the basis of nonlinear Recurrence Quantification 
Analysis 11 . It is important to emphasize that the selected channels are 
located either on the left-hemisphere or along the midline of the 
scalp. The midline sites in which we have observed hypnotizability 
related differences are closely linked to the default mode circuit, 
which includes the posterior cingulate cortex and precuneus, 
medial prefrontal and pregenual cingulate cortices, temporo- 
parietal regions, and medial temporal lobes (for a review see Hoeft 
et al. 24 ). This circuit is implicated in episodic memory retrieval, self- 
reflection, mental imagery, and stream-of-consciousness processing, 
which are cognitive activities potentially associated with long-lasting 
relaxation. The crucial role of these brain regions justified the 
selection of the four channels CZ, TP7, P3 and POl used in our 
analysis. Therefore the particular choice of electrodes was 
motivated by a compromise between the significance of the results 
and the simplicity of the overall procedure. 

Regarding signal acquisition and preprocessing, we emphasize 
that standard electroencephalographic recordings were performed 
according to the International 10-20 System by Ag/AgCl electrodes 
embedded in an elastic cap. The signal was sampled at 1 kHz. A 
band-pass digital filter (0.5, 100 Hz) was applied to EEG signals to 
remove DC bias and non- significant high frequency contributions; a 
notch filter was used to reduce the 50 Hz power-line noise. The 
FASTER algorithm 25 was used to remove the ECG and EOG con- 
taminations. It is based on well-known Independent Component 
Analysis (ICA) and removes undesired components by comparing 
the autocorrelation of the ECG and EOG signals to the correspond- 
ing independent ones. 

The EEG time series were obtained after applying a basic filtering 
procedure, starting from the real raw signal. There are two reasons 
why we used minimal manipulation of data: to avoid nonlinear 
methods, like those performed to obtain absolute values and/or 
amplitudes which could produce artificial time series; and, as we 
are predominantly interested in pre-clinical studies, our main task 
was to produce less elaborate and highly accessible procedures. 
Indeed, using a combined application of RQA and DFA techniques 
to our filtered time series, we obtained a useful biomarker in order to 
classify and distinguish lows from highs patients, which was the 
principal aim of our work. 

Combining recurrence and fluctuation analysis. In order to obtain 
a powerful discriminating procedure, we applied the DFA and RQA 
methods to EEG recordings relative to channels CZ, TP7, P3 and 
POl. The analysis was performed by using non-overlapping time 
windows of 10 seconds, which implies having analyzed subsets of 
10000 points of the main signal. For each channel and each 
individual we obtained a time series of fluctuation exponents a 
with a sampling rate of 10 seconds. We then set up a combined 
procedure to discriminate between EEG signals of highs and lows 
based on fluctuation exponent a and determinism D. In this way 
we obtained a dataset for each individual consisting of values of 
two measures calculated every 10 seconds for each of the four 
aforementioned channels. 

Fig. 1 and Fig. 2 report the resulting dataset for D and a, respect- 
ively, and for each of the channels CZ, TP7, P3 and POl. Each graph 
has been obtained through the evaluation of the average value of D 
and i indicators for highs and lows. In order to obtain a clearer 
discrimination of the two classes we also considered a smoothed 
(moving average) version of measures. On the left-hand sides of 
the same Figures we also included some examples of Recurrence 
Plot and log-log plots for the estimate of the indicators D and a. 
Moreover, we show in Table 1 mean and standard deviation of both 
indicators evaluated on channel CZ. 

The overall range of values of a runs from 1.1 to 1.3, proving 
nontrivial autocorrelation properties in EEG signals. In fact, the 
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Figure 1 | Examples of Recurrence Plots (left) obtained on samples of the EEG time series of two subjects involved in the experiments (see section 
Methods). Plots of Determinism extracted from EEG times series and corresponding to the most significant channels CZ, TP7, P3 and POl (right). 




Figure 2 | Examples of log-log plots for the estimate of the Fluctuation Exponent a (left). Plots of a extracted from EEGs and corresponding to the most 
significant channels CZ, TP7, P3 and POl (right). 
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range of values of the scaling exponent is crucial to the attempt to 
characterize EEG signals through their correlation properties. It is 
expected that these values depend on the state of participants, 
whether they are resting or involved in some cognitive tasks. 
Further, the same values depend on the kind of statistical analysis 
that is performed on the data. Our values are partly in agreement 
with other approaches. See Watters 26 , Linkenkaer-Hansen et al. 27 and 
Buiatti et al. 2 " for similar results. 

All values of D are high, showing significant levels of determinism 
in EEG signals. In general, D ranges from 0 to 1, and accounts for the 
distribution of the diagonal lines composed by recurrent points in the 
RP (see Figure 1 ) with respect to the total number of recurrent points. 
Noisy signals have D close to 0, while periodic signals have D close to 
1. Intermediate values may indicate time series produced by more 
complex systems, such as chaotic ones. 

Once the estimates for a exponents are added to the information 
extracted from measures of the deterministic part of the same signals 
we obtain complementary descriptions of the dataset concerning 
both the deterministic and the stochastic components. The couples 
of values (D, a) taken together provide a discriminating tool to dis- 
tinguish different levels of hypnotic susceptibility. 

Neural networks. The above results show a strong cross -evidence of 
a discrimination between highs and lows in the first seven minutes 
of analysis. In order to confirm this discrimination we employed 



Neural Networks, a very commonly used method for classification 
problems 29 . In this paper, the Neural Networks have been used to 
classify highs and lows on the basis of nonlinear measures, such as 
fluctuation exponent and determinism described above. For this 
purpose a two-layer feed-forward network, with sigmoidal hidden 
and output neurons, was used employing the Levenberg-Marquardt 
algorithm. We created a dataset of 576 vectors, each consisting of 8 
features (4 values of a, and 4 values of D) relating to each individual 
low and high, examined in a 10-second window during an overall 6- 
minutes time span, excluding the first transient minute. Thus the 
overall number of vectors is obtained from 36 temporal records times 
16 subjects. We used 70% of the whole dataset for training the 
network, 15% for validating and 15% for testing. We trained each 
network 1000 times, on both smoothed and non-smoothed data, for 
a number of neurons equal to 5, 10, 15, 20, respectively. All the 
networks have been trained on dataset related to each 10-second 
window, and then, having fixed the number of neurons and the 
type of dataset (smoothed or non-smoothed), we computed the 
average performance of the network as a function of time. 

The results obtained through these neural networks training are 
shown in Fig. 3, including both smoothed and non-smoothed signals. 
The time range used for the setup of the training dataset, and con- 
sequently the time interval shown in the graph, has been limited to a 
window of six minutes starting from the end of the first minute until 
the end of the seventh minute of our analysis. 

In Fig 4 the performances of the neural networks are highlighted 
through the evaluation of errors distribution. As one can see, the 
neural network with 10 neurons is the best choice, indeed the errors 
are significantly lower than those obtained by the neural network 
with 5 neurons. Furthermore, in both cases, the networks evaluated 
on smoothed data show better results. 

In order to additionally validate the outcomes of our experiment 
we included in Fig. 3 the results obtained from the neural network 
procedure applied to each indicator separately. It is clear that the 
combined use of the two indicators demonstrated much better per- 
formance compared to the neural network procedure based on each 
indicator separately. 
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Figure 3 | Time course of the errors obtained by different kinds of Neural Networks: 5 neurons (light blue), 10 neurons (red), 15 neurons (dark blue), 
20 neurons (green). The results for non-smoothed (dashed lines) and smoothed (continuous lines) signals, respectively, are highlighted. 
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Figure 4 | Histograms of the errors for the best (10 neurons) and worst (5 neurons) neural networks. 



Reliability and significance of results. Additional information on 
the method's performance can be extracted from Fig. 3, reporting the 
time course of the average errors for each group of neural networks. 
We find that for all groups the errors in the first minutes of the 
recordings are higher. Furthermore, initially the errors decrease 
and reach a flat region. This region was observed for all threshold 
settings and can be identified in the time interval ranging from 
second 170 to second 250 (see the green strip of Fig. 5). After this 
the errors increase again until the final part of the experiment (see 
also Fig. 3). 

Table 2 provides additional details of the above analysis, where the 
performances of the groups only focusing on the flat region ( 1 70-250 
seconds) are calculated. Specifically, the threshold in the first column 
and the number of selected networks falling below each threshold in 
the second column are reported. The third and fourth columns indi- 
cate the averages and the standard deviations of classification errors 
on the whole dataset and on the time interval 170-250 seconds, 
respectively. 

These results indicate that the identified neural networks show 
robustness with respect to the percentage errors. Indeed, the number 
of neural networks belonging to each group increases by raising the 
threshold and the same happens for the averages and standard devia- 
tions. Moreover, the number of networks for each group increases in 
a nonlinear way with respect to the threshold. The 25% of the total 
number of networks shows an error lower than 2.5% (first row of 
Table 2), while 90% of networks has an error lower than 20% (second 
last row of Table 2). A sort of saturating mechanism of networks 
number is also evident for threshold between 10% and 20%. Note also 
that performances improve significantly if we consider only the data 
in the time window 170-250 seconds. 

Discussion 

In this paper the problem of assessing hypnotic susceptibility has 
been approached using nonlinear dynamics through the combined 
use of Recurrence Quantification Analysis and Detrended 



Fluctuation Analysis techniques, supplemented by a neural network 
classification. By using the time course of two significant indicators, 
cross-evidence for hypnotic susceptibility was obtained starting from 
the basal EEG monitored on subjects in their ordinary state of 
consciousness. 

The main outcome of this paper is twofold: on one hand we 
obtained a clear indication on the EEG minimum recording time 
required to classify subjects as highs or lows in the absence of hyp- 
notic induction; on the other hand, the neural network procedure 
turned out to provide successful and robust classification of subjects, 
letting us take into account the intrinsic nonlinearities present in the 
signals. These results can be used by operators to design suitable 
experimental sessions related to hypnotic assessment, and imple- 
ment user-friendly devices, allowing quick hypnotic assessment in 
both experimental and clinical contexts. 

Thus, the combined use of the two techniques, RQA and DFA, 
inherited from nonlinear dynamics is a successful discriminating 
procedure and a hypnotic susceptibility classification tool, providing 
reliable predictions of hypnotizability levels. The strength of our 
approach relies on the use of two complementary features - one 
related to the deterministic part, and the other to the stochastic part 
of the signals - enabling to capture two different aspects of the com- 
plex nonlinear underlying dynamics. To the best of our knowledge, 
the integration of the two methods, although employed in different 
fields 3031 , has not been sufficiently used in EEG analysis. Some recent 
papers which make combined use of different classification methods 
based on different biomarkers measured from EEG signals are pre- 
sent in the literature regarding early diagnosis of Alzheimer's disease 
(AD). In particular, in Lehmann et al. 32 the authors used linear and 
non-linear classification algorithms - but not DFA - to discriminate 
between the electroencephalograms of patients with varying degrees 
of AD and their age-matched control subjects. Moreover in Poil 
et aV 3 the authors combine DFA together with other biomarkers 
in order to obtain a global diagnostic classification index in predict- 
ing the transition from normal aging to Alzheimer's disease. Both 
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Figure 5 | Time course of selected groups of neural networks' errors. 

approaches are crucial for effective application of early diagnostics 
and subsequent treatment strategies, but in neither of the previous 
examples the authors apply their classification methods combining 
RQA and DFA as indicators, nor their aim was to classify hypnotic 
susceptibility. In this respect our paper is quite innovative in the 
framework of hypnotic classification. 

With respect to the cited previous work carried out by some of the 
authors 11 , whose target was to perform a statistical analysis in order 
to obtain an initial understanding of possible discrimination between 
highs and lows subjects, the present paper is not a mere validation of 
previous results. In fact, in that paper discriminant analysis was 
applied after Bonferroni corrected MANOVA, showing significant 
differences between the highs and lows' determinism for the channels 
analysed in the present study. On the contrary, our main objective 
was to develop a clustering tool based on extended nonlinear indi- 
cators and on neural network models, which is reliable and robust 
enough to be used for clinical purposes. "Clustering tool" means that 



classification is automatic and does not require any additional ana- 
lysis. Once the tool is ready to use it becomes more and more reliable 
any time the EEC of a new subject is recorded and included in the 
system, due to the increase of the training and validating data sets of 
the neural network. 

Methods 

In this section the nonlinear methods of Detrended Fluctuation Analysis, Recurrence 
Quantification Analysis and Neural Networks are presented. 

Detrended Fluctuation Analysis. Detrended Fluctuation Analysis (DFA) can be 
traced back to the work by Hurst 34 , the so-called Rescaled Range (R/S) Analysis, 
introduced with the aim of quantifying fluctuations in time series, in order to 
distinguish random from fractal time series and to recognize the existence of long- 
range correlations. 

Roughly speaking, the R/S analysis gives a measure of displacement that the system 
undergoes - on average and rescaled by the local standard deviation - over the 
considered time. In particular, the procedure gives an estimate of exponent H in the 
scaling law for the mean quadratic displacement ((Ax) 2 ) ~ (At) 2H , by splitting a time 
series into adjacent non-overlapping windows and evaluating the range R of the 
standardized cumulative series rescaled by the standard deviation S. This scaling law 
is a generalization of the characteristic behavior in a Brownian motion, corresponding 
to H = 0.5, and it is related to fractional Brownian motion and anomalous diffusion. 

The original procedure by Hurst has been widely improved and generalized. Its 
most used extension is DFA, originally introduced to address the existence of power- 
laws in DNA sequences 35 . Through the specification of a proper detrending operation 
DFA proved itself to be a well suited procedure to obtain a correct estimation of 
power-law scaling, and a method specifically adapted to the analysis of non stationary 
time series 17,36 . 

To illustrate the DFA method, we consider a noisy time series, {x(i)} i= i N . We 
integrate the original time series and obtain the profile of cumulative sums {y(j)}: 

m= E wo -<*>]. <*>=ix>(0 f 1 ) 

N 

We divide the profile into N s — — boxes of equal size s. Since the length of the series is 

generally not an integer multiple of time scale s, a small portion of data at the end of 
the profile can not included in any interval. In order not to neglect this data we repeat 
the entire procedure starting from the end of the series, obtaining 2N S segments. Then, 
in each box we fit the integrated time series by using a discretized polynomial function 



Table 2 | Performances of neural networks for different levels of 
percentage errors (smoothed data). Column 1 reports the threshold 
used to identify the groups, column 2 the number of neural networks 
belonging to each group, column 3 the average and standard 
deviations of percentage errors on the total dataset and column 4 
the average and standard deviations of percentage errors in time 
interval 1 70-250 seconds 
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yfit, the local trend. We calculate the local trend for each interval of width 5 via a linear 
fit of data. Indicating the fit in the v-th box with^ f (i) we define the detrended 
profile as follows: 

^(i)=r(i) (v-l)s<i<vs (2) 

For each of 2N S intervals we calculate the mean square deviation from the local trend 



Yth-l)s+i] 



(3) 



Finally, we calculate the mean on all segments to obtain the fluctuation function 



F{s) 



A 



2iV s 



(4) 



The above computation is repeated for each box size (different scales) to provide a 
relationship between F(s) and s. The existence of a power-law like F(s) oc s a can be 
checked through the measure of the parameter a, called the scaling exponent, or 
correlation exponent - directly related to the Hurst Exponent H if the time series is 
stationary. It gives a measure of the correlation properties of the signal: if a — 0.5 the 
signal is an uncorrelated signal (white noise); if a < 0.5, the signal is anticorrelated; if a 
> 0.5, there are positive correlations in the signal 34,36 ' 37 . 

Recurrence Quantification Analysis. In this section we provide basic notions on 
recurrence methods for nonlinear time series data (for a review the reader is referred 
to Marwan et al. 13 ). 

Starting from a time series ...,s N ] where sj = s(iAt) and At is the sampling time, 
the system dynamics can be reconstructed using the Takens' embedding theorem 16 . 
The reconstructed trajectory x is expressed as a matrix in which each row is a phase 
space vector Xj — [s (r s J + T , . . . ,s j+ ^ Dli _i^] where i — 1, N — {D E — 1)t. The matrix 
is characterized by two key parameters: the embedding dimension D E and the delay 
time t. D E is the minimum dimension at which the reconstructed attractor can be 
considered completely unfolded and there is no overlapping of the reconstructed 
trajectories. The delay time t represents a measure of correlation existing between two 
consecutive components of D E - dimensional vectors used in trajectory 
reconstruction 38 . 

The Recurrence Plot (RP) is defined on the embedded trajectory x as follows: 



(5) 



where i and j are the time instant labels of the associated phase-space vector x. The RP 
is a two-dimensional object accounting for recurrences between the state vectors at 
different times. The recurrences are represented by black dots in the RP. Since any 
point is recurrent with itself, the RP always includes the diagonal line, R/j = 1, Vf — j, 
called Line of Identity (LOI). Thus the RP provides information on the dynamical 
properties of the reconstructed state space, with emphasis on the presence of peri- 
odicities of any length. 

Recurrence Quantification Analysis (RQA) provides a set of measures that allow to 
quantify the characteristics of RPs, otherwise difficult to be read by only visual 
inspection. In particular, RQA offers a set of indicators computed on the structures of 
the RP, based on the length / of the diagonal lines parallel to the LOI and their 
distribution P(l). This distribution provides valuable information regarding the 
structure of the RP and the unknown dynamics of the system under investigation 13 . 
Among the available indicators, Determinism D evaluates the proportion of recurrent 
points forming diagonal line structures in the RP and affects time series 
predictability 39 . 

In the present paper the RP for each electrode and each subject has been con- 
structed and the recurrence indicator Determinism has been quantified. In particular, 
the embedded state space has been reconstructed by using space dimension D E = 3 
and time lag t = 1200 40 , and D evaluated from its definition: 



(6) 



D accounts for the fraction of recurrent points forming diagonal structures with a 
minimum length l miri with respect to all recurrences. 

In order to provide a quantitative description of dynamical behavior of spatio- 
temporal systems starting from observed time series, such as the EEGs analyzed in 
this paper, RQA measures have been recently extended (see, for example Marwan 
et al. 41 , Facchini et al. 42 and Mocenni et al. 43 ). 

Neural Network Classification. Neural networks are computational models inspired 
by the nervous system of animals that are used in artificial intelligence applications, 
such as machine learning and pattern recognition. They are very commonly used for 
solving classification problems. 

Artificial neural networks provide procedures performing nonlinear statistical 
modeling, emulating biological neural systems. They give constructive alternatives to 
linear or logistic regression, the most commonly used methods for developing pre- 
dictive models for binary outcomes in diagnostic procedures. In this respect, neural 
networks can perform tasks that a linear program can not because they are able to 
catch the nonlinear features in the underlying system. They are simple and natural to 



setup, not very sensitive to noise and do not require complicated statistical training. 
Several training algorithms available in the literature can be used for their imple- 
mentation, in any kind of applications. The main disadvantage corresponds to their 
intrinsic "black box" nature, together with their basic empirical character, with the 
consequence that it is almost impossible to avoid the "human component" in clas- 
sification results. They do not always guarantee convergence towards an optimal 
solution. In our case the network was small and easy to construct. It worked very fast 
and well even if we recognize that the whole procedure could be implemented. Finally, 
we emphasize that the neural networks were only exploited in order to validate the 
hypothesis of distinction between highs and lows based on results obtained through 
completely different methodologies. 

Neural computing has recently emerged as a practical and successful technology, 
especially in problems of pattern recognition 44 . In this case it makes use of feed- 
forward network architectures such as multi-layer perceptron and the radial basis 
function network. Pattern recognition encompasses a wide range of information 
processing problems of great practical significance, in particular for medical dia- 
gnosis. In this respect neural networks are an extension of conventional techniques in 
statistical pattern recognition, and this field offers many powerful results. 

Basically simple artificial nodes ("neurons") are connected together to form a 
network which mimics a biological neural network. A neuron consists of a set of 
adaptive weights, which are numerical parameters tuned by a learning algorithm; they 
represent the strength of the connection between neurons activated during the 
training phase. At the end of the training the whole network is capable of approx- 
imating non-linear functions of their inputs 29 . 

Ethics statement. The work developed in this paper satisfies the ethics statement 
required by the journal. The experimental session was approved by the Ethics 
Committee of the University of Pisa where the EEG recordings were collected. 
Participants were selected from the students attending Physiology classes with 
informed written consent. 
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